PROBLEM 2.1 


Radiation is used to heat a plate of thickness L and thermal 
conductivity k. The radiation has the effect of volumetric energy 
generation at a variable rate given by 


q’'(x) = q7 e™ 


where q% and b are constant and x is the distance along the 


plate. The heated surface at x = O is maintained at uniform temperature T, while the 
opposite surface is insulated. Determine the temperature of the insulated surface. 


(1) Observations. (i) This is a steady state one-dimensional problem. (ii) Use a rectangular 
coordinate system. (iii) Energy generation is non-uniform. (iv) Two boundary conditions are 
needed. The temperature is specified at one boundary. The second boundary is insulated. (v) 
To determine the temperature of the insulated surface it is necessary to determine the 
temperature distribution in the plate. 


(2) Origin and Coordinates. A rectangular coordinate system is used with the origin at one 
of the two surfaces. The coordinate x is oriented as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) stationary and (4) constant 
properties. 


(ii) Governing Equations. Introducing the above assumptions into eq. (1.8) gives 


d°T q" 
+—=0 a 
D r (a) 


where 


k = thermal conductivity 

q" = rate of energy generation per unit volume 
T = temperature 

x= independent variable 


Energy generation q” varies with location according to 


q" =qze” (b) 


where q” and b are constant. Substituting (b) into (a) gives 
2 m 
T = 
aT „ooto (c) 
dx k 


(iii) Boundary Conditions. Since (c) is second order with a single independent variable 
x, two boundary conditions are needed. 


(1) Specified temperature at x = 0 
T(0) = T, (d) 
(2) The surface at x = L is insulated 
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PROBLEM 2.1 (continued) 


dT(L) _ 
dx 
(4) Solution. Separating variables and integrating (c) twice 


0 (e) 


q —bx 
T(x) =-— e "4+C,x4+C (f) 
bk 1 2 


where Cı and C2 are constants of integration. Application of boundary conditions (d) and 
(e) gives Cı and Cz 


do -bL do 
C =-= e”, C =T, +~ 
1 bk 2 o b?k (2) 
Substituting (g) into (f) 
T(x) = = [i-e -bxe™ |+7, (h) 


The temperature of the insulated surface is obtained by setting x = L in (h) 
T(t) =—[1-a+ bbe |+7, (i) 
b°k 
Expressing (h) and (i) in dimensionless form gives 


alae =O Dalle (j) 


and 


PAS he. =1-(1+bL)e (k) 


m 
oO 


bk 
(5) Checking. Dimensional check: Each term in solution (h) has units of temperature. 
Boundary conditions check: Setting x = 0 in (h) gives T(0) =T,. Thus boundary condition 


(d) is satisfied. To check condition (e) solution (h) is differentiated to obtain 


= fe [oe — be] 


(m) 


Evaluation (m) at x = L shows that condition (e) is satisfied. 


Differential equation check: Direct substitution of (h) into equation (c) shows that the 
solution satisfies the governing equation. 


Global energy conservation check: The total energy generated in the plate, £ g» must be 
equal to the energy conducted from the plate at x = 0. Energy generated is given by 
: L L m 
E, = fa" ax = af e™dx = l-e”) (k) 
0 0 


Energy conducted at x = 0, E(0), is determined by substituting (m) into Fourier’s law 
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PROBLEM 2.1 (continued) 


E(0) __,~- GO _ 
dx 


de. -bL 
-=(l-e `) () 
b 
The minus sign indicates that heat is conducted in the negative x-direction. Thus energy is 
conserved. 


Limiting check: For the special case of q” =0, temperature distribution should be uniform 
given by7T(x)=T,. Setting g” =0 in (h) gives the anticipated result. 


(6) Comments. (i) The solution to the special case of constant energy generation 
corresponds to b=0. However, setting b=0 in solution (h) gives 0/0. To obtain the 
solution to this case it is necessary to set b=O in differential equation (c) and solve the 
resulting equation. (11) The dimensionless solution is characterized by a single dimensionless 
parameter bL. 
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PROBLEM 2.2 


Repeat Example 2.3 with the outer plates generating energy at a rate q" and no energy is 
generated in the inner plate. 


(1) Observations. (i) The three plates form a composite wall. (ii) The geometry can be 
described by a rectangular coordinate system. (iii) Due to symmetry, no heat is conducted 
through the center plate. Thus the center plate is at a uniform temperature. (iv) Heat 
conduction is assumed to be in the direction normal to the three plates. (v) This problem 
reduces to a single heat generating plate in which one side is maintained at a specified 
temperature and the other side is insulated. 


(2) Origin and Coordinates. Since only one plate needs to be analyzed, the origin and 
coordinate x are selected as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state and (3) constant conductivity of 
the heat generating plate. 


(ii) Governing Equations. Let the subscript 2 refer to the heat generating plate. Based 
on the above assumptions, eq.(1.8) gives 


T, 
d T, m x 
2 —=0 (a) 
dx k, 
Where EN EEE IT 
k, = thermal conductivity 
q" = rate of energy generation per unit volume o T, 


T = temperature 
x = coordinate 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
x, two boundary conditions are needed. 


(1) Specified temperature at x = 0 


T,(0) =T, (b) 
(2) Insulated surface at x = L, 
dT, (Ly) _ 0 © 
dx 


where L, is the thickness of the heat generating plate. 


(4) Solution. Integrating (a) twice gives 


m 


T, (x) ss 


x’ +C x+C, (d) 
2 


where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives C, and C, 
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PROBLEM 2.2 (continued) 


"T, 
G=f2, C2 =T, (e) 


2 


Substituting (e) into (d) gives the temperature distribution in the heat generating plate 


m 2 


T(x) = 4 [2L /x)-1]+T, (f) 
2k, 


To determine the temperature distribution in the center plate we apply continuity of 
temperature at x = L, 


T,(L,) =T,(L,) (g) 


Noting that the center plate must be at a uniform temperature, evaluating (f) at x = L, and 


substituting into (g) gives 
my2 


L 
TW=N) =)= +T, (h) 
2 
Expressing (f) and (h) in dimensionless form gives 
T,(x)-T, 1x . 
m -= 2 (x / L, )] (i) 
ql, 2h, 
k, 
and 
T, (x) =i; Z 1 (j) 
gl 2 
ky 


(5) Checking. Dimensional check: Each term in (f) has units of temperature. 


Boundary conditions check: Substitution of solutions (f) into (b) and (c) shows that the two 
boundary conditions are satisfied. 


Differential equation check: Direct substitution of (f) into equation (a) shows that the 
solution satisfies the governing equation. 


Global energy conservation check: The total energy generated in the plate, E,, must be 


£ kj 
equal to the energy conducted from the plate at x = 0. Energy generated is given by 


È, =q"AL, (k) 


where A is surface area. Energy conducted at x = 0, E(0), is determined by differentiation 
(f) and substituting into Fourier’s law 


E(0) =-k, 


A dT, (0) = —q"AL, (m) 
dx 


The minus sign indicates that heat is conducted in the negative x-direction. Equations (k) 
and (m) show that energy is conserved. 
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PROBLEM 2.2 (continued) 


Limiting check: If energy generation vanishes, the composite wall should be at a uniform 
temperature T,. Setting q” =Oin (£) and (h) gives T,(x) =7,(x) =T,. 


(6) Comments. (i) An alternate approach to solving this problem is to consider half the 
center plate and the heat generating plate as a composite wall, note symmetry, write two 
differential equations and four boundary conditions and solve for T, (x) and T,(x). Clearly 
this is a longer approach to solving the problem. (ii) No parameter characterizes the solution. 
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PROBLEM 2.3 


A plate of thickness L, and conductivity k, moves with a T, 
velocity U over a stationary plate of thickness L, and i 
conductivity k,. The pressure between the two plates is P 
and the coefficient of friction is u. The surface of the 
stationary plate is insulated while that of the moving plate is 
maintained at constant temperature T,. Determine the steady 
state temperature distribution in the two plates. 


LASS ASS SSS SSS 


(1) Observations. (i) The two plates form a composite wall. (ii) The geometry can be 
described by a rectangular coordinate system. (iii) Heat conduction is assumed to be in the 
direction normal to the two plates. (iv) Since the stationary plate is insulated no heat is 
conducted through it. Thus it is at a uniform temperature. (v) All energy dissipated due to 
friction at the interface is conducted through the moving plate. (vi) This problem reduces to 
one-dimensional steady state conduction in a single plate with specified flux at one surface 
and a specified temperature at the opposite surface. (vii) Plate motion plays no role in the 
governing equation since no heat is conducted in the direction of motion. 


(2) Origin and Coordinates. Since only one plate needs to be analyzed, the origin is 
selected at the interface and the coordinate x is directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant conductivity of the 
moving plate, (4) no energy generation and (5) uniform energy dissipation at the interface. 


(ii) Governing Equations. Let the subscript 1 refer to the moving plate. Based on the 
above assumptions, eq.(1.8) applied to the moving plate gives 


=0 (a) 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
x, two boundary conditions are needed. 


(1) Specified heat flux at x = 0 . Fourier’s law gives 


dT, (0 
=i anw) _ MPU (b) 
dx 
where 
k, = thermal conductivity 
P = interface pressure 
U = plate velocity 
u= coefficient of friction 
(2) Specified temperature at x = L, 
T(L,) =T, (c) 


where L, is the thickness of the moving plate. 
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PROBLEM 2.3 (continued) 
(4) Solution. Integrating (a) twice gives 


where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives C, and C, 


PU 


Substituting (e) into (d) gives the temperature distribution in the moving plate 


uPUL į 


T(x) = 1-(x/L,)]+T, (f) 


1 


To determine the temperature distribution in the stationary plate we apply continuity of 
temperature at x = 0 


T,(0) =7T,() (g) 


Noting that the stationary plate must be at a uniform temperature, evaluating (f) at x = 0 and 
substituting into (g) gives 


HES if (h) 


T,(x) =T, (0) =7, (0) 7 o 
1 


Expressing (f) and (h) in dimensionless form gives 


eo =|—-— (i) 
MPUL, L 
kı 
T(x) -T, ; 
oO =] 
HMPUL, Q) 
kı 


(5) Checking. Dimensional check: Each term in (f) has units of temperature 


Boundary conditions check: Substitution of solutions (f) into (b) and (c) shows that the two 
boundary conditions are satisfied. 


Differential equation check: Direct substitution of (f) into equation (a) shows that the 
solution satisfies the governing equation. 


Limiting check: If no energy is added at the interface (u =0or U = 0 or P = 0), the two 
plates should be at a uniform temperature T,. Setting any of the quantities 4, U or P equal 
to zero in (f) and (h) gives 7,(x) =T,(x) =T,,. 


(6) Comments. (i) An alternate approach to solving this problem is to consider the two 
plates as a composite wall, write two differential equations and four boundary conditions 
and solve for T, (x) and T,(x). Clearly this is a longer approach to solving the problem. (ii) 
No parameter characterizes the solution. 
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PROBLEM 2.4 


A thin electric element is wedged between two plates of 
conductivities k, and k,. The element dissipates uniform 
heat flux q”. The thickness of one plate is L, and that of 
the other is L,. One plate is insulated while the other 
exchanges heat by convection. The ambient temperature is 
T,, and the heat transfer coefficient is h. Determine the 


temperature of the insulated surface for one-dimensional 
steady state conduction. 


(1) Observations. (i) The two plates form a composite wall. (ii) The geometry can be 
described by a rectangular coordinate system. (iii) Heat conduction is assumed to be in the 
direction normal to the two plates. (iv) Since no heat is conducted through the insulated 
plate, its temperature is uniform equal to the interface temperature. (v) All energy dissipated 
by the element at the interface is conducted through the upper plate. (vi) This problem 
reduces to one-dimensional steady state conduction in a single plate with specified flux at 
one surface and convection at the opposite surface. (vii) To determine the temperature of the 
insulated surface it is necessary to determine the temperature distribution in the upper plate. 


(2) Origin and Coordinates. Since only one plate needs to be analyzed, the origin is 
selected at the interface and coordinate x is directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant conductivity of the 
upper plate, (4) no energy generation and (5) uniform energy dissipation at the interface. 


(ii) Governing Equations. Let the subscript 1 refer to the upper plate. Based on the 
above assumptions, eq.(1.8) applied to the upper plate gives 
d°T, 
F = 
dx 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
x, two boundary conditions are needed. 


(a) 


(1) Specified heat flux at x = 0 


=k; ae =q; (b) 
(2) Convection at x = L,. Fourier’s law and Newton’s law of cooling give 
-y FU = hR) -T] Q 
where L, is the thickness of the upper plate. 
(4) Solution. Integrating (a) twice gives 
T(x) =C,x+C, (d) 
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PROBLEM 2.4 (continued) 


where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives C, and C, 


C=-2, C, = [hL /k)+1]+T, (e) 
1 


Substituting (e) into (d) and rearranging gives the temperature distribution in the upper plate 


_ qo | hL x 
T, (x)= a a Dular, (f) 


To determine the temperature distribution in the insulated plate we apply continuity of 
temperature at x = 0 


T (0) =7, (0) (2) 


Noting that the stationary lower plate must be at a uniform temperature, evaluating (f) at x = 
0 and substituting into (g) gives 


1 


n hL 
T, (x) = T, (0) = T, (0) -gH siler, (h) 


Equation (h) gives the temperature of the insulated surface. 


Expressing (f) and (h) in dimensionless form gives 


T, -T 
i) (i 
do Li 
h 
and 
T, -T 
O le 214 Bi (i) 
LA 
h 
where Bi is the Biot number defined as 
Bi = hh (k) 
kı 


(5) Checking. Dimensional check: Each term in (f) has units of temperature. 


Boundary conditions check: Substitution of solution (f) into (b) and (c) shows that the two 
boundary conditions are satisfied. 


Differential equation check: Direct substitution of (f) into equation (a) shows that the 
solution satisfies the governing equation. 


Limiting check: (i) If no energy is added at the interface, the two plates should be at a 
uniform temperature T. Setting q% = 0 in (f) and (h) gives T,(x) =T,(x) =T.. (ii) If h = 0, 
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PROBLEM 2.4 (continued) 


no heat can leave the plates and thus their temperature will be infinite. Setting h = 0 in (£) 
and (h) gives T,(x) = T, (x) =. 


Qualitative check: Increasing interface energy input increases the temperature level in the 
two plates. Equations (f) and (h) exhibit this behavior. 


(6) Comments. (i) An alternate approach to solving this problem is to consider the two 
plates as a composite wall, write two differential equations and four boundary conditions 
and solve for T (x) and T,(x). Clearly this is a longer approach to solving the problem. (ii) 
The problem is characterized by a single dimensionless parameter AL /k, which is the Biot 
number. This parameter appears in problems involving surface convection. 
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PROBLEM 2.5 
A thin electric element is sandwiched between two plates of h, 
conductivity k and thickness L each. The element dissipates a o 
flux qù. Each plate generates energy at a volumetric rate +L Ue gr I 
of q” and exchanges heat by convection with an ambient fluid L d 
at T,,. The heat transfer coefficient is h. Determine the temperature hT, 


of the electric element. 


(1) Observations. (i) The two plates form a composite wall. (ii) The geometry can be 
described by a rectangular coordinate system. (iii) Heat conduction is assumed to be in the 
direction normal to the two plates. (iv) Symmetry requires that half the energy dissipated by 
the element at the interface is conducted through the upper plate and the other half is 
conducted through the lower plate. Thus only one plate needs to be analyzed.(v) The 
problem reduces to one-dimensional steady state conduction in a single plate with energy 
generation having specified flux at one surface and convection at the opposite surface. (vi) 
To determine the temperature of the element it is necessary to determine the temperature 
distribution in one of the two plates. 


(2) Origin and Coordinates. The upper plate is considered for analysis. The origin is 
selected at the upper surface and coordinate x is directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) LT, 
constant conductivity of the upper plate, (4) uniform energy T m 
generation and (5) uniform energy dissipation at the interface. L X o — 

(ii) Governing Equations. Let the subscript 1 refer to the Dh 
upper plate. Based on the above assumptions, eq.(1.8) applied to 2 
the upper plate gives 

dT m 
144 9 (a) 
dx k 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
x, two boundary conditions are needed. 


(2) Specified heat flux at x = L 


T L n 
k d ne ) = do (b) 
dx 2 
(2) Convection at x =0. Fourier’s law and Newton’s law of cooling give 
T, 
EAE O] () 
dx 
where Lis the thickness of the upper plate. 
(4) Solution. Integrating (a) twice gives 
Hoe- P26 226 (d) 
2k 
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PROBLEM 2.5 (continued) 


where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives C, and C, 


"n "I, 
; C=. 42242 (e) 


“ K 2h h 


Substituting (e) into (d) and rearranging gives the temperature distribution in the upper plate 


n "T, m 
1.4 3 2 


q! q"L 
T, (x) =T, + + + x 
1(4) =T, E i T (f) 


2h h 


Expressing (f) in dimensionless form gives 


2 
T,(x)-T. i : 
pF) py Go ljepi | (g) 
q"Llh 2q"L L] 2|L 


where Biis the Biot number defined as 


hL 
Bi = — h 
5 (h) 
To determine the temperature of the electric element set x = L in (g) 
T,(L) -T. i j 
1E) -To __4o hepi] (i) 
q"LIh 2q"L 2 


(5) Checking. Dimensional check: Each term in (g) is dimensionless. 


Boundary conditions check: Substitution of solution (f) into (b) and (c) shows that the two 
boundary conditions are satisfied. 


Differential equation check: Direct substitution of (f) into equation (a) shows that the 
solution satisfies the governing equation. 


Limiting check: (i) If no energy is added at the interface and there is no energy generation 
the two plates should be at a uniform temperature T. Setting q% =q” =0 in (f) gives 
T,(x) =T,,. (ii) If h = 0, no heat can leave the plate and thus the temperature will be infinite. 
Setting h = 0 in (f) gives T,(x) =. 


Qualitative check: Increasing interface energy input and/or energy generation increases the 
temperature level in the plate. Equation (f) exhibit this behavior. 


(6) Comments. (i) An alternate approach to solving this problem is to consider the two 
plates as a composite wall, write two differential equations and four boundary conditions 
and solve for T,(x) and T,(x). Clearly this is a longer approach to solving the problem. (ii) 
The solution is characterized by two parameters: the convection parameter Bi and the ratio of 
energy dissipated in the element to energy generation, q” /q"L. 


0 
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PROBLEM 2.6 


One side of a plate is heated with uniform flux q” while the other side exchanges heat by 
convection and radiation. Surface emissivity is €, heat transfer coefficient h, ambient 
temperature T,,, surroundings temperature T „, plate thickness L and the conductivity is k. 


Assume one-dimensional steady state conduction and use a simplified radiation model. 
Determine: 


EE es ped gir 
[a] The temperature distribution in the plate. 


hT,, 
[b] The two surface temperatures for: 
L X k 
h=27W/m?—°C, k=15W/m-°C, L=8cm, € = 0.95, t 
q! =19,500Wim?, T.,,= 18°C, T,,= 22 °C. tT FOF gh 


(1) Observations. (i) This is a one-dimensional steady state problem. (ii) Two boundary 
conditions are needed. (iii) Cartesian geometry. (iv) Specified flux at one surface and 
convection and radiation at the opposite surface. (v) Since radiation is involved in this 
problem, temperature units should be expressed in kelvin 


(2) Origin and Coordinates. The origin and coordinate x are shown. 
(3) Formulation. 


(i) Assumptions. (1) Steady state, (2) one-dimensional, (3) constant properties, (4) no 
energy generation, (5) small gray surface which is enclosed by a much larger surface, (6) 
uniform conditions at the two surfaces and (7) stationary material. 


(ii) Governing Equations. Based on the above assumptions, the heat equation in 
Cartesian coordinates is given by 


aT 
dx? 
(iii) Boundary Conditions. Two boundary conditions are needed. They are: 


(1) Specified flux at x =0 


0 (a) 


dT) _, 


-k 
dx do 


(b) 


where 
k =thermal conductivity = 15 W/m—°C 
q! = surface flux = 19,500 W/m? 


(2) Convection and radiation at x=L. Conservation of energy, Fourier’s law of 
conduction, Newton’s law of cooling and Stefan-Boltzmann radiation law give 


-ZO D-T +e oit o-ri,] i 
x 
where 


h =heat transfer coefficient = 27 W/m? —° C = 27 W/m? -K 
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PROBLEM 2.6 (continued) 

L = plate thickness = 8 cm 

T „= ambient temperature = 22 °C = (22 + 273.15) K = 295.15 K 

T.,,. = Surroundings temperature = 18°C = (18 + 273.15) K = 291.15 K 

g = emissivity = 0.95 

o = Stefan-Boltzmann constant = 5.67 x108 W/m? —K4 
(4) Solution. Integrating (a) twice gives 

T(x) =C,x4+C, (d) 

where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives 


do 
g =i-Ér+c; -1, reo] C Ér+c -Ti | (f) 


Equation (f) gives the constant C,. However, this equation must be solved for C, by trial 
and error. 


(5) Checking. Dimensional Check: Each term in equations (b), (c) and (f) has units of flux. 


Differential equation check: Solution (d) satisfies the governing equation (a). 


(6) Computations. Equation (e) gives 


2 o 
G= 19,500 (W/m) _ 1300 C _ 300% 
15 (W/m-°C) m m 
Equation (f) gives 


- 19,500 (W/m?) 
15 (W/m-K) 


19,500(W/m?) = 27 (W/m? - 9 0.08(m) +C, — 2951500) 


-19,500 (W/m°) 
15 (W/m -K) 


4 
+ (0.95)5.67 x10 (W/m? - «| 0.08(m) + cs} z esns) 
Solving this equation for C, by trial and error gives 
C, =761.9K 
Substituting C, and C, into (d) gives 
T(x) =761.9(K)—1300(K/m)x (g) 


Surface temperatures at x = 0 and x = L are determined using (g) 
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PROBLEM 2.6 (continued) 
T(0) =761.9K 
T(L) = 657.9 K 
(7) Comments. (i) A trial and error procedure was required to solve this problem. (ii) Care 
should be exercised in selecting the correct units for temperature in radiation problems. (ili) 


With the two temperatures determined a numerical check based on conservation of energy 
can be made. 


Energy flux entering plate = Energy flux conducted through the plate = Energy flux leaving (h) 
Energy flux entering plate = q” =19,500 W/m? 


T(O)-T(L) 


Energy flux conducted through the plate = Lk 


Energy flux leaving = i|-Eisc, -1, [+eo|(-#L+¢,)4 Ti, 


Numerical computations show that each term in (h) is equal to 19,500 W/m’. 
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PROBLEM 2.7 


A hollow shaft of outer radius R, and inner radius R; rotates h,T, 
inside a sleeve of inner radius R, and outer radius R,. Frictional 
heat is generated at the interface at a flux q%. At the inner shaft R 
surface heat is added at a flux qf. The sleeve is cooled by 
convection with a heat transfer coefficient h. The ambient 
temperature is T,,. Determine the steady state one-dimensional 
temperature distribution in the sleeve. qn 


sleeve 
(1) Observations. (i) The shaft and sleeve form a composite 


cylinder. (ii) Use a cylindrical coordinate system. (i11) Heat conduction is assumed to be in 
the radial direction only. (iv) At steady state, energy entering the shaft at R; must be equal 
to energy leaving at R,,. (v) Frictional energy dissipated at the interface and energy added at 
R; must be conducted through the sleeve. (vi) This problem reduces to one-dimensional 
steady state conduction in a hollow cylinder (sleeve) with specified flux at the inside surface 
and convection at the outside surface. 


(2) Origin and Coordinates. The origin is selected at the center and the coordinate r is 
directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant conductivity of the 
sleeve, (4) no energy generation, and (5) uniform frictional energy at the interface. 


(ii) Governing Equations. Based on the above assumptions, eq.(1.11) applied to the sleeve 
gives 
—(r—) =0 (a) 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
r, two boundary conditions are needed. 


(1) Specified heat flux at r = R, . Fourier’s law gives 

dT(R , 
s) = ds (b) 
dr 


where q% is the heat flux added to the sleeve at r=R,. This flux is the sum of frictional 
energy generated at the interface and energy added at R;. That is 


-k 


qa =4; +q, (c) 
where 


m 


q) = heat flux leaving the shaft at R, 
q; = heat flux generated by friction at the interface R, 
q", = total heat flux entering the sleeve at R, 


At steady state, conservation of energy applied to the shaft requires that energy entering the 
shaft at R, be equal to energy leaving it at R,. Thus 
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PROBLEM 2.7 (continued) 


2nRiq) = 20R,q, 


R. 
"= iq" d 
do n“ (d) 


S 


Substituting (c) and (d) into (b) gives the boundary condition at R, 


dT(R) „ R, 
=q 


—k +— q? e 
i stR Gi (e) 
(2) Convection at r = R, . Fourier’s law and Newton’s law of cooling give 
dT(R,) 
-k TE = AlP(R,)-T] @) 


(4) Solution. Integrating (a) twice gives 
T(r) =C, Inr+C, (g) 
where Cand C, are constants of integration. Application of boundary conditions (e) and (f) 


gives C, and C, 


R, R, 
Ci =- las +R IR aT], Co =T + Lat + (RR ai K/AR,) +1nR,] (hy) 
Substituting (h) into (g) and rearranging gives the temperature distribution in the sleeve 
R 
P(r) =a; + (R; /R, ai Illa, 17) + (KIAR, )]+ T, (i) 
Expressing (1) in dimensionless form gives 


T-T q'R, | 


+—— |!In(R,/r)+(1/Bi j 
GR, q'R, [ ( o ) ( )] 0) 
k 
where the Biot number Bi is defined as 
hR 
Bi = — k 
Fi (k) 


(5) Checking. Dimensional check: Each term in (i) has units of temperature. Each term in (j) 
is dimensionless. 


Boundary conditions check: Substitution of solutions (i) into (e) and (f) shows that the two 
boundary conditions are satisfied. 


Differential equation check: Direct substitution of (i) into (a) shows that the solution 
satisfies the governing equation. 
Limiting check: (i) If no energy is added to the shaft at R; and no frictional heat is generated 


at the interface R,, the sleeve should be at a uniform temperature T,,. Setting q; =q% =0 in 
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PROBLEM 2.7 (continued) 


(i) gives T(r) =T.,. Gi) If h = 0, no heat can leave the sleeve and the temperature will be 
infinite. Setting h = 0 in (i) gives T(r) =. 

(6) Comments. (i) An alternate approach to solving this problem is to consider the shaft and 
sleeve as a composite cylinder, write two differential equations and four boundary 
conditions and solve for the temperature distribution in both shaft and cylinder. Clearly this 
is a longer approach to solving the problem. (ii) Two parameters characterize the solution: 
The Biot number and heat ratio (q/R;/q!R,). 
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PROBLEM 2.8 
A hollow cylinder of inner radius r, and outer radius r, is heated io 
with flux qat its inner surface. The outside surface exchanges fr h,T, E, 


heat with the ambient and surroundings by convection and 
radiation. The heat transfer coefficient is h, ambient temperature 
T„, surroundings temperature T,,, and surface emissivity is €. 
Assume steady state one-dimensional conduction and use a 
simplified radiation model. Determine: 


[a] The temperature distribution in the cylinder. 
[b] The inner and outer surface temperatures for: 


q! = 35,500Wim?, T,,= 24°C, T,,= 14°C, k=3.8W/m-°C, r, = 12cm, 


SUr 


r =5.5cm, h=31.4Wim’-°C, € = 0.92 


(1) Observations. (i) This is a one-dimensional steady state problem. (ii) Two boundary 
conditions are needed. (iii) Cylindrical geometry. (iv) Specified flux at the inner surface and 
convection and radiation at the outer surface. (v) Since radiation is involved in this problem, 
temperature units should be expressed in kelvin 


(2) Origin and Coordinates. The origin and coordinate r are shown. 
(3) Formulation. 


(i) Assumptions. (1) Steady state, (2) one-dimensional, (3) constant properties, (4) no 
energy generation, (5) small gray surface which is enclosed by a much larger surface, (6) 
uniform conditions at the two surfaces and (7) stationary material. 


(ii) Governing Equations. Based on the above assumptions, the heat equation in 
Cartesian coordinates is given by 


d ; dT 
dr a (a) 
(iii) Boundary Conditions. Two boundary conditions are needed. They are: 
(2) Specified flux at r = r, 
T(r; 
= k d (7;) = 7" (b) 


dr = 
where 

k =thermal conductivity = 3.8 W/m—°C 

q! = surface flux = 35,500 W/m? 


r; =inside radius = 5.5 cm 


(2) Convection and radiation at r=r,. Conservation of energy, Fourier’s law of 
conduction, Newton’s law of cooling and Stefan-Boltzmann radiation law give 
dT 
Ae) <afry,)-T.}+e0,)-Ta,] ©) 
F 


C 


-k 
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PROBLEM 2.8 (continued) 

where 

h = heat transfer coefficient = 31.4 W/m?—°C = 27 W/m? - K 

r, =outside radius = 12 cm 

T,,= ambient temperature = 24 °C = (24 + 273.15) K = 297.15 K 

T „„ = surroundings temperature = 14°C = (14 + 273.15) K = 287.15 K 

g = emissivity = 0.92 

o = Stefan-Boltzmann constant = 5.67 x108 W/m? — K4 
(4) Solution. Integrating (a) twice gives 

T(r) =C,Inr+C, (d) 

where C,and C, are constants of integration. Application of boundary conditions (b) and 
(c) gives 


qi 
c =-fn (©) 
gn) Lr in, +C, -1, reo] (Lin, +C,)* -rá (f) 


Equation (f) gives the constant C,. However, this equation must be solved for C, by trial 
and error. 


(5) Checking. Dimensional Check: Each term in equations (b), (c) and (f) has units of flux. 


Differential equation check: Solution (d) satisfies the governing equation (a). 


(6) Computations. Equation (e) gives 


2 o 
35,500 (W/m ) 0.055) = 513.82 -= 513.82 É 
m 


3.8 (W/m-° C) m 


1 — 
Equation (f) gives 


0.055(m) - 0.055 (m) 


35,500 (W/m?) =31.4 (W/m? - K| 35,500(W/(W/m? ) In 0.12(m).+ Cy — 297.15 «| 


0.12(m) 3.8 (W/m - K) 
-0.055(m) i 
+(0.92)5.67 x108 (Wim? - xt fesse 35,500 (W/m ) In 0.12(m) + c,| (287.15) 


Solving this equation for C, by trial and error gives 
C, =-484.8K 
Substituting (e) and C, into (d) gives 


T(r) = 484.8 (K) - ay Inr (g) 
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PROBLEM 2.8 (continued) 


To determine the inside surface temperature set r = r, in (g) 


35,500(W /m?—°C) 
3.8(W /m—°C) 


T(r,) =—484.8(K) 0.055(m) 1n0.055(m) = 1005.48K 


Similarly, the outside surface temperature is 


35,500(W /m?—°C) 
3.8(W/m—°C) 


T(r, ) = —484.8(K) 0.055(m) In0.055(m) = 604.625K 


(7) Comments. (i) A trial and error procedure was required to solve this problem. (ii) Care 
should be exercised in selecting the correct units for temperature in radiation problems. (ili) 
With the two temperatures determined a numerical check based on conservation of energy 
for a tube of length L can be made: 


Energy entering tube at r; = Energy conducted through tube = Energy leaving tube at r, (h) 
Energy entering tube = - =2a7r,q' Wim 


T(r,)-T 
Energy conducted through tube = ; = (i= TU) 


1 r, 
— jn 
2rk b 


l 


Energy flux leaving at r, = a =27 r hT, )-T,|+27r,e olric) -T4 | 


SUr 


Numerical computations show that each term in (h) is equal to 12267.9 W/m. 
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PROBLEM 2.9 


Radiation is used to heat a hollow sphere of inner radiusr,, outer radius r, and 


conductivity k. Due to the thermal absorption characteristics of the material the radiation 


results in a variable energy generation rate given by 
2 
5 
q” (r) = q” = 
oO 


where q” is constant and r is the radial coordinate. The 


inside surface is insulated and the outside surface is 
maintained at temperature T,. Assume steady state one- 
dimensional conduction, determine the temperature 
distribution. radiation 


(1) Observations. (i) Radiation results in a non-uniform 

energy generation. (ii) Use a spherical coordinate system. (iii) At steady state the total 
energy generated in the sphere must be equal to the energy leaving the surface by convection 
and radiation. (iv) This is a one-dimensional steady state conduction problem in a hollow 
sphere with insulated inner surface and convection and radiation at the outer surface. 


(2) Origin and Coordinates. The origin is selected at the center and the coordinate r is 
directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant conductivity, (4) 
uniform surface and ambient conditions, and (5) the sphere is modeled as a small gray 
surface which is enclosed by a much larger surface. 


(ii) Governing Equations. Based on the above assumptions, eq.(1.13) gives 


1 dy .dTy q" 
r +—=0 a 
r? m a k ( ) 
where q” depends on r according to 
2 
m m r 
q"r)=q > (b) 
ro 
Substituting (b) into (a) 
l1 dọ dT q" r’ 
r + =0 c 
r? I = k r? (© 


(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
r, two boundary conditions are needed. 


(1) Insulated boundary at 7; 


aT(7;) -0 (d) 
dr 
PROBLEM 2.9 (continued) 
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(2) Specified temperature at r, 
T(r,)=T, (e) 


(4) Solution. Separating variables and integrating (a) twice gives 


r- Į_î 


r +C (£) 
20kr2 r” 
where C, and C, are constants of integration. Application of boundary condition (d) gives 
n 5 
do ti 
C ==> 
153k (g) 
Boundary condition (2) gives 
n2 5 
doro |l ti 
C, =T, + -+= (h) 
oe Sk E | 
Substituting (g) and (h) into (f) 
" 5 5 
T(r) =T, -2 rt +4% -ri -44 (i) 
20kr; r o 
Expressing (i) in dimensionless form gives 
T-T, 1 fn tlet a olel ð 
q'r? 20 5 i Sin & 20| 


k 


(5) Checking. Dimensional check: Each term in (i) has units of temperature and each term 
in (j) is dimensionless. 


Boundary conditions check: Equation (i) satisfies boundary conditions (d) and (e). 


Differential equation check: Direct substitution of (i) into (c) shows that the solution 
satisfies the governing equation. 


Limiting check: If no energy is generated in the sphere the temperature will be uniform equal 
toT,. Setting g? =0 in (i) gives T(r) =T,,. 


oO 
Conservation of energy check: The total energy generated in the sphere must leave the outer 
surface by conduction. Total energy generate is obtained by integrating (b) over the volume 
of the hollow sphere. Energy conducted the outer surface is obtained by applying Fourier’s 
law at r = r, using solution (i). Both are found to be equal to 
An do, 5 


Total energy generated = g2 Co r) 
r 


oO 


(7) Comments. According to (j), the solution is characterized by a single dimensionless 
geometric parameter r; / r; . 
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PROBLEM 2.10 


An electric wire of radius 2 mm and conductivity 398 W/m—°C generates energy at a rate 
of 1.25x10° W/m}. The surroundings and ambient temperatures are 78°C. and 82°C, 
respectively. The heat transfer coefficient is 6.5 W/m?—°C and surface emissivity is 0.9. 
Neglecting axial conduction and assuming steady state, determine surface and centerline 
temperatures. 


(1) Observations. (i) The wire generates energy. (ii) Use a 

cylindrical coordinate system. (iii) Heat conduction is assumed to T 
be in the radial direction only. (iv) At steady state the total energy Pa 
generated in the wire must be equal to the energy leaving the 

surface by convection and radiation. (v) To determine center and E 
surface temperatures it is necessary to solve for the temperature 

distribution in the wire. (vi) This is a one-dimensional steady state 

conduction problem in a solid cylinder with energy generation. s 


S~ 
N 


(2) Origin and Coordinates. The origin is selected at the center 
and the coordinate r is directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant conductivity, (4) 
uniform energy generation, (5) uniform surface and ambient conditions, and (6) the wire is 
modeled as a small gray surface which is enclosed by a much larger surface. 


(ii) Governing Equations. Based on the above assumptions, eq.(1.11) gives 


ld, dT q" 
=0 
mr a j k (a) 


where 

k = thermal conductivity = 398 W/m—°C 

q" = energy generation rate = 1.25 x10° W/m? 

(iii) Boundary Conditions. Since (a) is second order with a single independent variable 
r, two boundary conditions are needed. 

(2) Symmetry gives 

dT(O) _ 
dr 


0 (b) 


(2) Convection and radiation at the surface r= r,. Conservation of energy at the surface 
and using Fourier’s law, Newton’s law of cooling and Stefan-Boltzmann radiation law give 


- 10) n6p¢,)-1,]ee0[T*(,)-Tey] E 
r 


where 


h =heat transfer coefficient = 6.5 W/m? —°C 
r, = wire radius = 2 mm = 0.002 m 
T „„ = surroundings temperature = 82°C = 355 K 
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PROBLEM 2.10 (continued) 


T,, =ambient temperature = 78°C = 351 K 
e =emissivity = 0.9 
o = Stefan-Boltzmann constant =5.67x 10° W/m? — K* 
(4) Solution. Separating variables and integrating (a) twice gives 


m 


q 
To)=-7p" +Cilnr+C, (d) 
where C, and C, are constants of integration. Application of boundary condition (b) gives 
C =0 (e) 
Thus (d) becomes 
q” 
T(r)=C,- = (f) 
Boundary condition (2) gives 
Penk le, -(@"r2 14k) - 7, |+ eo |C, -(q"r2140)} -c 0oT4, (g) 


This equation gives C,. However, it can not be solved explicitly for C,. Solution must be 
obtained by trial and error. Once C, is determined, equation (f) gives center and surface 
temperatures by setting r=0 and r = r,, respectively. Thus 


T(0)=C, (h) 
and 
T(r,)=C,-Lr (i) 
4k 


(5) Checking. Dimensional check: Each term in (g) has units of flux and each term in (i) has 
units of temperature. 


Boundary conditions check: Equation (f) satisfies boundary condition (b). Condition (c) can 
not be checked since the constant C, is not known explicitly. 


Differential equation check: Direct substitution of (f) into (a) shows that the solution 
satisfies the governing equation. 


Limiting check: If no energy is generated in the wire and T, =T,,,, the wire temperature 


will be uniform equal to the ambient temperature. Setting g” =0 and T, =T,,, in (f) gives 
T(r) =C,. Equation (g) gives C, = T}. 


(6) Computations. Substituting numerical values in (g) and solving for C, by trial and error 
gives C, =361.2K. Equations (h) and (i) give center and surface temperatures as 


T(0)=361.2K and 7(r,)=361.2K 
(7) Comments. The wire is essentially at uniform temperature. This is due to the fact that k 


is high and r, and h are small. 
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PROBLEM 2.11 


The cross-sectional area of a fin is A,(x) and its circumference is C(x). Along its upper half 
surface heat is exchanged by convection with an ambient fluid atT,,. The heat transfer 
coefficient is h. Along its lower half surface heat is exchanged by radiation with a 
surroundings atT,,,.. Using a simplified radiation model formulate the steady state heat 


equation for this fin. 


(1) Observations. (i) This is a fin 
problem with surface convection and 
radiation. (ii) To formulate the fin 0 
equation apply conservation of energy 

to an element of the fin using fin 
approximation and Fourier, Newton 
and Stefan-Boltzmann laws. 


(2) Origin and Coordinates. The 
origin is selected at one end of the fin 
and the coordinate x is directed as 
shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) no energy generation, (4) 
Bi <<1, and (5) the fin is modeled as a small gray surface which is enclosed by a much 
larger surface. 


(ii) Governing Equations. To formulate the heat equation for this fin we select an 
element dx at location x as shown. This element is enlarged showing its geometric 
characteristics and energy interchange at its surfaces. Conservation of energy for the element 
under steady state conditions requires that 


Ee = E out (a) 
where 


E,, = rate of energy added to element 


in 


È jy, = rate of energy removed from element 


Energy enters the element by conduction at a rate q, and leaves at arate q,+ (dq, /dx) dx. 
Energy also leaves by convection at a rate dq,and by radiation at a rate dg,. We now 
formulate each term in (a) 


Ein= 4, (b) 
and 
: d qx 
F out = Ix + —~ dx +.dq, + dq, (c) 
dx 
Substituting (b) and (c) into (a) 
LT +dq,+dq, =0 (d) 
dx 
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PROBLEM 2.11 (continued) 


We introduce Fourier, Newton and Stefan-Boltzmann laws to eliminate q,, dq, and dq,, 
respectively. Thus, 


=—-kA = (e) 
qd» a Cc dx e 
dq, =h(T -T „)dA,/2 (f) 
and 
dq, =co|r* 13, A, 12 (2) 


where A,=A.(x) is the cross-sectional area through which heat is conducted. The 
infinitesimal area dA, /2 is half the surface area of the element where heat is exchanged by 
convection at the top half and by radiation at the bottom half. 


Substituting (e), (f) and (g) into (d) we obtain 


A, 
= (ka, w2 2- h(T - T) -E0 (Tt -Tar) 5 Z (h) 
The gradient of surface area dA, /dx is given by eq. (2.6a) 
= = C(x)[l+ (dy, dx)? }” (i) 
X 


where y,(x) is variable which describes the fin profile and C(x) is the circumference. If k 
is assumed constant, (h) becomes 


2 
dT 1ldAdT_ _ h dA, (r-T,)-£2 dA, (rt -T4.)=0 D 
dx? A, dx dx 2kA, dx 2kA, dx 


(4) Checking. Dimensional check: Each term in (j) has the same units. 


Limiting check: If there is no radiation, (j) should reduce to eq. (2.5b). Setting £ =0 in (j) 
gives eq. (2.9). 


(5) Comments. (i) Radiation introduces non-linearity in the fin equation. (ii) To obtain the 
corresponding equation for a fin with no surface convection set h = 0 in (j) and remove the 
2 factor if radiation takes place over the entire surface. 


https://gioumeh.com/produ@2$heat—conduction-solution/ 


PROBLEM 2.12 


A constant area fin of length 2L and cross-sectional area A, generates heat at a volumetric 


m 


rate q”. Half the fin is insulated while the other half exchanges heat by convection. The 
heat transfer coefficient is h and the L < L 

ambient temperature is T,,. The base and j j h, To 

tip are insulated. Determine the steady of x q" ; 
state temperature at the mid-section. At LT 


2- 00 


what location is the temperature highest? 


(1) Observations. (i) This is a constant area stationary fin. (ii) Temperature distribution can 
be assumed one dimensional. (iii) Surface conditions are different for each half. Thus two 
equations are needed. (iv) Both ends are insulated. (v) The highest temperature is at the 
insulated end x = 0. 


(2) Origin and Coordinates. The origin is selected at one end of the fin and the coordinate 
x is directed as shown. 


(3) Formulation. 


(i) Assumptions. (1) One-dimensional, (2) steady state, (3) constant cross-sectional area, 
(4) uniform energy generation, (5) Bi <<1, (6) constant conductivity, (7) no radiation and 
(8) uniform h and T.,. 


(ii) Governing Equations. Since half the fin is insulated while the other half exchanges 
heat by convection, two fin equations are needed. Let 


T, (x)= temperature distribution in the first half, O< x < L 
T, (x)= temperature distribution in the second half, L<x<2L 


The fin equation for each section is obtained from eq. (2.5b) 


2 1A m 
es a a (2.5b) 
dx” A.(x) dx dx kA,.(x) dx k 


Note that A.is constant and dA, /dx = C , where C is the circumference. Eq. (2.5b) becomes 


d 2T hC q" 

me a, ( ) ; (a) 
In the first half the fin is insulated along its surface and thus the convection term vanishes. 
Equation (a) gives 


d°T, q” 
+——=0 O<x<L b 
dx? k 9 ( ) 
For the second half (a) gives 
dT. m 
= (T, Ta) + Č =0, L<x<2L (c) 
where 
2 AC 
m“ =— d 
a (d) 
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PROBLEM 2.12 (continued) 


(iii) Boundary Conditions. Four boundary conditions are needed. 


(1) Insulated surface at x = 0 


ORN 5 
dx 
(2) Equality of temperature at x = L 
TO =T,(L) (f) 


(3) Equality of flux at x = L 


dT (L) _ dT,(L) 


dx dx (8) 
(4) Insulated surface at x = 2L 
dT (2L) o i 
dx 
(4) Solution. Separating variables and integrating (b) twice 

q” 2 , 

T, (x) =-— x" +Ax+B i 
1x) 7) (i) 


where A and B are constants of integration. The solution to (c) is given in Appendix A, 
equation (A-2b) 


T(x) =T,, + A > +Csinhmx + Dcoshmx (j) 
m 


where C and D are constants of integration. Application of the four boundary conditions 
gives the four constants 

A=0 (k) 

B-T + ge N q”  q"L tanh 2mLsinhmL — coshmL 

“2k km? km tanh 2mLcoshmL —sinhmL 


D 


q"L tanh 2ml 
C= - (m) 
km sinhmL — tanh 2mLcoshmL 


jati — (n) 
km tanh 2mLcoshmL — sinhmL 


Substituting into (i) and (j) gives the solutions as 


m2 m m . m 
T(x) =T, + q E E: : q”L tanh 2mLsinhmL cone q 2 (0) 
2k km km tanh2mLcoshmL—sinhmL 2k 


and 
T,(x) =T,, + q : +4 L oe ae 
km km tanh 2mLcoshmL — sinh mL 


(p) 


https://gioumeh.com/produés3#heat—conduction-solution/ 


